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<N ■ Abstract 

Spatial structure of Cooper pairs with quantum numbers color 3*, I = J = 
L = S = in ud 2 flavor quark matter is studied by solving the gap equation 
and calculating the coherence length in full momentum range without the 
J> ' weak coupling approximation. Although the gap at the Fermi surface and 

the coherence length depend on density weakly the shape of the r-space pair 
ly-j ■ wave function varies strongly with density. This result indicates that quark 

Cooper pairs become more bosonic at higher densities. 
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Color superconducting phases of strongly interacting matter at high density and low 
temperature []l],|2||| are attracting much attention recently. They were studied first as an 
example of pair condensation in relativistic many-body systems by Bailin and Love in the 
early 80's |4j] (see also Ref. ||). They also mentioned the nucleon-nucleon pairing as an- 
other example; its detailed study was begun in the early 90's || and is developing recently 
|7U8l]^,|T0[1 . A study of the quark-quark pairing in an SU(2) color model was done at almost 
the same time [O (see also Refs. fC^JIBfl ). Iwasaki and Iwado's work on a 1 flavor system |TJJ 
was the first study of the realistic SU(3) color system (see also Ref. |H|). Since the works of 
Rapp et al. [IE] and Alford et al. |T7| on ud 2 flavor system, color superconductivity has been 
studied extensively. The 1 5'o ( J — L — S — 0) state in the 2 flavor case |18|] and the color- 
flavor- locked state in uds 3 flavor case Hl9| ,p0i^lip2i p3l , respectively, have been understood 



to be the most favored channels. The magnetic interaction has been shown to be responsible 
for these pair condensations |2^j25|j26|j27|j28|j29|] . Their astrophysical consequences were also 
studied ^^. 

The purpose of this Brief Report is to visualize the spatial structure, in particular its 
density dependence, of quark Cooper pairs, which has not been discussed, to the author's 
knowledge. This is done for Cooper pairs of the simplest form; color 3*, isosinglet 1 Sq pairs 
in the 2 flavor system. To this end, calculations are performed for zero temperature under an 
instantaneous approximation while full k-dependence is retained. That is, we allow strong 
coupling in the sense that momenta far from the Fermi surface also contribute. We will 
mention briefly the frequency dependence later. 

We start from deriving the gap equation for A(k) representing the gap for color 3*, 
isosinglet 1 So pairs composed of two quarks which are time- reversal conjugate to each other 
with respect to space and spin, 



Ak s /i,kV/v — (— I) 2 s 5k-kiS s - s ieffieu>A(k) 



(i) 



where s, / and i denote spin, flavor and color, respectively. Note that time-reversal is repre- 
sented by T times complex conjugate, and T = —iC'y 5 . We work in a Hamiltonian formalism 



since it is more convenient than the Gor'kov formalism |32| in the present study in which 
only the pair condensate is considered. The latter might be more suitable when the coupling 



between the pair condensate and the chiral condensate |L7| is considered. Such a formalism 
has already been developed for the nuclear system 0Q] although the N-N condensate is 
negligibly small M due to a large nucleon mass. The resulting 12 x 12 Hamiltonian matrix, 



(E k -n A \ 

^ -A* -{E k -n)J ' 

E k = ^k 2 + M|, f i = E kF , 

is easily block-diagonalized to two 6x6 matrices by inspection 
agonalized along the lines of the 1 flavor case studied in Ref 
condensate in terms of the Bogoliubov coefficients, we obtain 



(2) 

Then they are fully di- 
ll]. By expressing the pair 
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As for the residual interaction, we adopt the one gluon exchange interaction with the 
leading order screening. Although the gluon propagator therein contains a gauge-dependent 
term in general, it vanishes due to the equation of motion of the external quark |27j . Then 



we proceed to a static (instantaneous) approximation |I3, I5|jI8| . At this step, the dynamic 
magnetic screening drops. By performing a spin average and an angle integration to project 
out the S'-wave component, we obtain 
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The running coupling constant is taken from Ref. |33 as in Refs. [15,18 
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with parameters q 2 = 1.5Aq CD , Aq CD = 0.4 GeV. 

The frequency dependence of the gap, which is beyond the present formulation, has 
been shown recently to change the asymptotic behavior of the gap in the weak-coupling 
approaches p4| , p6| , p8| . There frequency dependence appears both in the gluon propagator 
that determines the interaction and in the other part of the integrand in the Eliashberg 
equation p^ , |26| , p3| . Only the effects of the former can be estimated approximately here by 
modifying the magnetic part of Eq. (ffl) as 
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where m 2 ^ is originating from the Landau damping p6| , p7| and the frequency q° is an in- 
dependent parameter here due to the nature of the present formulation. Evidently this 
simply weakens the magnetic attraction and accordingly monotonically reduces the gap at 
the Fermi surface, A(g° = 10 MeV)/A(g° = 0) = 0.70, A(g° = 50 MeV)/A(g° = 0) = 
0.53, at k-p = 2 fm _1 , for example, while the g°-dependence of the coherence length (defined 



later) is rather weak. In the Eliashberg equation the gap is given as a function of p°, where 
q° = p° — k° and the equation includes the integration with respect to k°. Its p°-dependence 
was numerically studied in Refs. |26],|23| , where the gap exhibits a maximum at intermediate 
densities whereas decreases monotonically at asymptotically high densities. 

In the following, we concentrate on the static case and present numerical results in three 



steps. As for the quark mass, M q = 10 MeV is adopted according to Ref. [|Tj]. First we 
discuss the gap at the Fermi surface A(kp) and the coherence length £ as functions of the 
Fermi momentum kp, which is related to the baryon density as p = 2k F /3ir 2 in the present 
2 flavor case. The coherence length is defined as [33|, 

1/2 

(6) 



Jo \dk\ K aK 



in the strong coupling case, in terms of the pair wave function, 

which is identical to the pair condensate up to phase factors. In Fig.|l|(a), A(kp) is graphed 
for kp = 1.5 - 3.25 fm -1 , which corresponds to p/po — 1-5 - 15 if the normal density of 
symmetric (JV = Z) nuclear matter is defined as po = 2(kp)l/3ir 2 , (&f)o = 1-30 fm -1 |35| . 
This shows very weak fcp-dependence; the superconducting phase survives up to practically 
infinite density; for example, A(kp = 60 fm -1 ) ~ 50 MeV. We confirmed that the gap 
is predominantly brought about by the magnetic interaction; about 50% of the total gap 
remains when the electric interaction is cut artificially while only 5 - 10% remains when 
the magnetic part is cut. This magnetic dominance can also be deduced from the strong 
/cp-dependence of v(kp, k) (Fig.f2](c) below) except at k ~ k F where the magnetic interaction 
gives very strong attraction irrespective of kp, in contrast to the weak /cp-dependence of 
A(fep). The coherence length is shown in Fig.|l|(b). This can be compared with the Pippard 
length in the weak coupling theory, £ = kp/nA(kp)p ~ 1/7tA(A;f) (see for example, Ref. 
pfifl); the magnitudes of £ are approximated fairly well by £o whereas the /cF-dependence is 

1 l"\ 

different. The average interparticle distance d = \\) /^f, derived from p q = 3p = l/d 3 , 
is also shown in the figure. The magnitudes of £ and d are very similar to each other as in 
the nucleon-nucleon case |1(J , whereas £ is 3 - 4 orders of magnitude larger than d in metals. 



This indicates strong coupling feature in the sense that quark Cooper pairs are compact and 
therefore bosonic. That is, fermion exchange is less when their mutual overlap is small. If 
M q is changed to 100 MeV and 300 MeV, for example, A(/c F ) decreases to 90 - 98 % and 
54 - 81 % of the original values, respectively, for k-p = 1.5 - 3.25 fm _1 , whereas the changes 
in C, are negligibly small. 

Another quantity which can be compared with £ is the London penetration depth Al, 



which is defined by J p/Ane^ps in the relativistic case, where e s and p s are the electric charge 
of the Cooper pair and the density of the superconducting component, respectively. This 
expression can be derived from the relativistic Ginzburg-Landau theory developed in Ref. 
ID (see also Refs. |30|j25[|). In the present ud 2 flavor case, Al — 20 - 8 fm for k F = 1.5 



- 3.25 fm 1 is obtained by adopting e s = |e + (— |)e and p s = |p ? . Although Bailin and 
Love concluded that quark matter is a superconductor of the first kind, £ > V2\l, based 



on their estimate, A(fcp) ~ 1 MeV [|J, the present quantitative study shows that quark 



matter is a superconductor of the second kind at these densities as pointed out in Ref. |30 
At asymptotically high densities, however, it changes to a superconductor of the first kind 
since the penetration depth decreases rapidly as density increases. One should note that 
these discussion applies only to a part of magnetic field since the original photon and a 
gluon combine to a massless 'rotated' photon although both the electromagnetic U(l) and 
the color SU(3) break in color superconductors as recently discussed in Ref. |3"7| ]. 

Next, we turn to fc-dependence at each kp. The compactness of the Cooper pairs men- 
tioned above indicates spreading of A(/c) (Fig.^|(a)) and <p(k) (Figj^(b)) in /c-space. Figure 
0(a) shows that A(k) spreads to larger k at higher densities as expected, while the width 
of 4>(k) in Fig.|2|(b), which corresponds to l/£, is almost /^-independent as mentioned above 
since the quasiparticle energy E'(k) in the denominator in Eq.(^) grows as k goes away from 
kp. The asymmetric shapes of <fi{k) reflect the smallness of £. Its width is similar to that of 
v(k F , k) in Fig.^J(c). 

Finally, we Fourier-transform <p{k) to 

<P{k) Jo {kr)k 2 dk, (8) 
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where jo{kr) is a spherical Bessel function, in order to look into the spatial structure of 
quark Cooper pairs more closely. This quantity has already been discussed for nucleon 
Cooper pairs in non-relativistic |]39| , |40[| and relativistic |10| studies. Although the coherence 



length £, corresponding to the root mean square distance with respect to this wave function, 
is almost /^-independent, the shape of <p(r) is strongly /cp-dependent; large-fc components 
in high-density cases bring nodes in r-space. Conversely, <p(r) is diffuse at low densities. In 
other words, quark Cooper pairs become more bosonic at higher densities. 

To summarize, we have studied numerically the spatial structure of quark Cooper pairs 
in the color 3*, isosinglet 1 So channel by solving the gap equation in full momentum range. 
Although the long-range magnetic interaction is predominantly responsible for the quark- 
quark pairing, the gap function spreads also in momentum space. The resulting coherence 
length is almost density independent as well as the gap at the Fermi surface and is of magni- 
tudes similar to the average interquark distance. The dependence of the pair wave function 
(condensate) on the relative momentum and on the distance between two quarks that form 
a Cooper pair has been presented. This indicates that quark Cooper pairs become more 
bosonic at higher densities although the coherence length is almost density independent. 



The author thanks T. Hatsuda for informing him of R. Horie's master thesis in Ref. [TE| 
and for the communication of Ref. [^ 
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FIG. 1. (a) Pairing gap at the Fermi surface, and (b) coherence length and average interquark 
distance, as functions of the Fermi momentum kp. 
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FIG. 2. (a) Gap function, (b) pair wave function, and (c) matrix element v(kp, k), as functions 
of the momentum k, calculated at kp = 1.5 - 3.0 fm _1 . Note that, in (c), k = k-p is excluded since 
at this point the diverging magnetic interaction does not contribute to the gap equation 
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FIG. 3. Pair wave function as a function of the distance r, calculated at kp = 1.5 - 3.0 fm 
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